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Abstract 


Bovine coronavirus isolates associated with recent outbreaks of respiratory disease in Ontario and Quebec dairy 
farms were compared to reference strains known to be responsible for neonatal calf diarrhea (NCD) or winter 
dysentery (WD) of adult cattle. In respect to their hemagglutinating properties and their higher RDE activities with 
rat erythrocytes, WDBCoV strains differed from NCDBCoV strains and respiratory bovine coronaviruses RBCoV 
strains. Serologically, three MAbs directed to the HE glycoprotein of the WDBCoV strain BCQ.2590 recognized two 
serogroups amongst NCDBCoV strains by hemagglutination inhibition, whereas only one of the MAbs failed to react 
toward three of the four RBCoV isolates tested. Sequencing analysis of the S (S1 portion), HE, ORF4 and ORF5 
genes of BCoV isolates associated with different clinical syndromes indicated that neither insertions or deletions could 
explain their distinct tropism. For the HE glycoprotein, a total of 15 amino acids (aa) substitutions were identified 
by comparing field isolates to the prototype Mebus strain. Two specific proline substitutions were identified for 
virulent strains being located in the signal peptides (aa 5) and aa position 367; one specific aa change was revealed 
at position 66 for RBCoV field isolates. Analysis of the S1 portion of the S glycoprotein revealed a total of eight aa 
changes specific to enteropathogenic (EBCoV) strains and eight aa changes specific to RBCoV strains. For all BCoV 
isolates studied, the region located between the S and M genes (ORF4) apparently encodes for two non-structural (ns) 
proteins of 4.9 and 4.8 kDa. A specific non-sense mutation was identified for the nucleotide at position 88 of the 
putative 4.9 kDa protein gene of RBCoV isolates resulting in 29 rather that 43 aa residues. The ORF5, which encodes 
a 12.7 ns protein and the 9.5 kDa E protein, was highly conserved amongst the BCoV field isolates. © 2001 Elsevier 
Science B.V. All rights reserved. 
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Bovine coronavirus (BCoV) is a member of the 
family Coronaviridae, order Nidovirales (De 
Vries et al., 1997). The genome of coronaviruses 
consists in a single stranded, capped, and 
polyadenylated positive-sense RNA molecule of 
approximately 26 to 30 kb in length, that yields a 
3’-coterminal nested set of subgenomic mRNAs in 
the infected cells (Spaan et al., 1981; Lai, 1986). 
Coronaviruses are mostly spherical in shape, en- 
veloped with a diameter ranging between 120 and 
150 nm, and possess an helical nucleocapsid. The 
virion contains four structural proteins: the spike 
(S) glycoprotein of 160-200 kDa which forms the 
long club-shaped surface projections, the matrix 
glycoprotein (M) of 20—38 kDa, the phosphory- 
lated nucleocapsid protein (N) of 50-60 kDa and 
the small membrane protein (E) of about 9.5 kDa 
(Spaan et al., 1988 (a review)). Based on serologi- 
cal cross-reactivities and genomic relatedness, the 
coronaviruses have been separated into three dis- 
tinct subgroups (Schockley et al., 1987; Dea et al., 
1990). The BCoV belongs to the second subgroup, 
with hemagglutinating encephalomyelitis virus of 
swine (HEV), human respiratory coronavirus 
(HCoV-0C43), turkey enteric coronavirus 
(TCoV) and murine hepatitis coronaviruses 
(MHVs). Members of this subgroup possess an 
additional 140 kDa envelope-associated glyco- 
protein, the hemagglutinin—esterase (HE), corre- 
sponding to a second fringe of granular surface 
projections (King et al., 1985; Hogue and Brian, 
1986; Dea et al., 1990). Both the S and HE 
glycoproteins are able to hemagglutinate red 
blood cells by binding to N-acetyl-9-O neuraminic 
acid (Vlasak et al., 1988; Schultze et al., 1991b), 
elicit the production of virus neutralizing antibod- 
ies (Dea et al., 1990; Vautherot et al., 1992), and 
are probably involved in determining the viral 
tissue and host tropism (Spaan et al., 1988). The 
HE is further associated to acetyl esterase (AE) 
and receptor destroying enzyme (RDE) activities 
(Vlasak et al., 1988; Parker et al., 1990). Few 
investigators have also postulated that the regions 
situated between the S and M genes (ORF4 and 
ORFS), potentially encoding non-structural (ns) 
proteins of 4.9, 4.8 and 12.7 kDa, and the small 
9.5 kDa membrane protein (E), could be involved 
in the presumed preferential respiratory tropism 


of HCoV-OC43 and HEV (Abraham et al., 
1990b; Mounir and Talbot, 1993; Vieler et al., 
1996). 

Althought BCoV isolates are generally recog- 
nized as being associated with enteric diseases of 
newborn calves (NCD: neonatal calf diarrhea) 
(Mebus et al., 1973) and winter dysentery (WD) 
or chronic diarrhea in adult cattle (Crouch et al., 
1985; Benfield and Saif, 1990), there exist much 
evidences that other isolates may be responsible 
or contribute to acute respiratory problems in 
feedlot cattle (Reynolds et al., 1985; Storz et al., 
1992). A recent survey conducted in feedlots or 
livestock shows in several states in the USA, 
emphasized the role of RBCoV in the bovine 
shipping fever syndrome pneumonia (Storz et al., 
1996, 2000). Despite their association to different 
clinical pathologies, BCoV isolates apparently be- 
long to a single serotype on the basis of virus 
neutralisation (VN) and hemagglutination inhibi- 
tion (HAI) tests using polyclonal hyperimmune 
sera. However, data obtained from recent serolog- 
ical findings using anti-HE and anti-S MAbs sug- 
gest the existence of different subgroups of BCoV, 
each containing EBCoV and WDBCoV strains 
responsible for either NCD, WD or chronic shed- 
ding in adult cattle (El Ghorr et al., 1989; Tsune- 
mitsu and Saif, 1995; Milane et al., 1997; 
Fukutomi et al., 1999). Antigenic and genomic 
variations have been also recently reported among 
EBCoV and RBCoV strains (Chouljenko et al., 
1998; Hasoksuz et al., 1999). The significance of 
such variability in the virulence and tropism of 
field isolates of BCoV is still unknown. 

The purpose of the present study was to analyse 
the biological, serological and genomic properties 
of cytopathogenic RBCoV strains recently iso- 
lated in Quebec and Ontario, Canada, in compari- 
son to the reference Mebus strain. We further 
compared the RBCoV strains to previously char- 
acterized Quebec EBCoV strains isolated in new- 
born calves with diarrhea and to Quebec 
WDBCoV strains in order to identify specific 
strain markers that should be considered for diag- 
nosis and development of vaccines. 

The cell culture-adapted and attenuated Mebus 
strain of BCoV was obtained from the American 
Type Culture Collection (ATCC VR_ 874, 
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Rockville, MD). This prototype NCD strain was 
originally isolated from diarrhea fluid of a 2- 
week-old calf and attenuated following more than 
30 successive passages in fetal bovine kidney 
(FBK) cells (Mebus et al., 1973). The 67N strain 
of porcine hemagglutinating encephalomyelitis 
virus (VR 741) was also obtained from the 
ATCC, whereas the HCoV-OC43 virus was 
kindly provided to us by PJ Talbot, INRS-Institut 
Armand-Frappier, Laval, Quebec, Canada. The 
RBCoV reference strain OK-0514, isolated from 
the nasal swab of a calf with clinical respiratory 
tract disease subsequent to shipping from Okla- 
homa to a feedlot in Kansas, was kindly provided 
to us by J. Storz, School of Veterinary Medicine 
of Louisiana (Storz et al., 1996). Field strains of 
EBCoV associated with either diarrhea in new- 
born calves or WD in adult cattle were obtained 
from affected dairy herds from different geo- 
graphical areas in Quebec, Canada (Dea et al., 
1995; Milane et al., 1997). The RBCoV strain 
BCQ.3994 was isolated from the lungs of pneu- 
monic growing calves from a dairy farm in South- 
ern Quebec that developed signs of respiratory 
tract disease during 1998, whereas two other 
Canadian RBCoV_ strains (BCO.43277 and 
BCO.44175) were isolated from pneumonic grow- 
ing calves with pneumonia from two different 
Ontario dairy farms during 1997 (kindly provided 
to us by S. Carman, Ontario Ministry of Agricul- 
ture and Fisheries, Ontario Veterinary College, 
Guelph, Ontario, Canada). All field isolates of 
BCoV were propagated for less than five passages 
in the G clone of humal rectal tumor cells (HRT- 
18G) developed through selection and medium 
modulation (Storz et al., 1996). Trypsin-treatment 
of field isolates, as well as routine addition of 10 
U/ml of bovine pancreatic trypsin in the mainte- 
nance medium of infected cultures, were done as 
described elsewhere (Dea et al., 1990). All the 
BCoV isolates propagated in the HRT-G cells 
behaved as highly fusogenic (induced production 
of small syncytia with less than 10 nuclei) and 
cytolytic (complete destruction of the monolayers 
within 4—5 days after infection) strains (data not 
shown). 

Hemagglutination (HA) of rat erythrocytes at 4 
and 37°C, and hemagglutination inhibition (HAI) 


tests were performed, as previously described 
(Dea et al., 1980). The receptor destroying (RDE) 
activities of the various BCoV isolates were deter- 
mined according to Storz et al. (1992), and ex- 
pressed as the ratio between HA titers of infected 
supernatant fluids determined after incubation at 
1 h at 4°C and HA titers determined after shifting 
to 37°C for 1 h, leading to activation of the 
enzyme activity and elution of the viral particles 
from the surface of the erythrocytes (Vlasak et al., 
1988). Following three freeze and thaw cycles of 
the infected cultures, the extracellular virions of 
the various BCoV isolates were purified by isopyc- 
nic ultracentrifugation on sucrose gradients (Dea 
et al., 1989). To determine the esterase activities 
associated to the HE glycoprotein of the various 
BCoV isolates, aliquots (15 pl) of purified viral 
preparations, adjusted to contain approximately 
the same amounts of viral particles as determined 
by negative stained electron microscopy (Dea et 
al., 1989), were added to 1 ml of PBS containing 
1 mM p-nitrophenyl acetate (PNPA). Hydrolysis 
of the substrate following an incubation of 5 min 
at room T° was monitored by determining the 
optical density at a wavelength of 405 nm (A4o5) 
(Parker et al., 1990). 

For serological studies, monoclonal antibodies 
1D6-2, 1D6-3 and 9F2-1R directed to the HE 
glycoprotein of the Quebec WDBCoV reference 
strain BCQ.2590 were used (Milane et al., 1997), 
as well as rabbit hyperimmune sera produced to 
the concentrated and purified Mebus and 
BCQ.2590 strains (Dea et al., 1995). The anti-HE 
MAbs were obtained following fusion experi- 
ments with splenocytes from Balb/c mice which 
have been immunized with an enriched peplom- 
eric fraction of the WDBCoV BCQ.2590 strain 
(Milane et al., 1997). The specificities of the 
MAbs toward 65 and/or 140 kDa polypeptides, 
corresponding to the monomeric and dimeric 
forms of the HE glycoproteins, have been deter- 
mined by Western blotting and radiotmmunopre- 
cipitation experiments. No reactivity was 
observed toward the S (100 or 200 kDa) glyco- 
protein, and none of these anti-HE MAbs had 
neutralizing activity against the homologous and 
Mebus BCoV strains (Milane et al., 1997). 
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The capacity of BCoV field isolates to attach to 
9-O-acetylated sialic residues present at the sur- 
face of rat erythrocytes at different incubation T° 
is a biological property which has been previously 
reported for distinguishing between NCD- and 
WD.-associated BCoV strains (Dea et al., 1995; 
Milane et al., 1997). In agreement with previous 
findings, HA titers of all BCoV strains studied at 
4°C varied between 640 and 5120, and only HA 
titers of WDBCoV strains diminished significantly 
when switching the incubation T° from 4 to 37°C, 
due to the activation of the higher RDE activity 
of the WDBCoV strains (Table 1). For RBCoV 
and NCDBCoV strains, the capacity to destroy 
the host cell receptors was not influenced by 
variation of the incubation 7°, except for 
BCQ.1523 strain for which a minimal RDE activ- 
ity of 2 was determined, meaning that its HA 
titers at 4°C was at least twice the value obtained 
at 37°C. The esterase activities of the WDBCoV 
and RBCoV strains were comparable and ap- 
peared superior to that of NCD strains. Study of 
the cross-reactivity of the anti-HE MAbs revealed 
a lack of reactivity (HAI titers <20) of MAbs 
1D6-2, 1D6-3 and 9F2-1R toward NCD-associ- 
ated Mebus and BCQ.571 strains, whereas HAI 
titers of at least 2560 were obtained when testing 
MAbs 1D6-2 and 1D6-3 against all other EBCoV 
and RBCoV strains studied. However, MAb 9F2- 
1R reacted only weakly (HAI titers of 40-80) to 
both Ontario and Quebec RBCoV strains tested, 
a titer of 2560 being obtained when tested against 
the reference OK-0514 RBCoV strain. Despite the 
lack of reactivity to the anti-HE MAbs tested, all 
BCoV strains tested reacted positively by HAI 
against both rabbit hyperimmune sera (anti- 
Mebus and anti-WD.2590). 

For the preparation of genomic RNA, aliquots 
(50 ul) of purified virus or infected cell lysates 
were supplemented with 1 ul of RNA guard 
(Pharmacia) and extraction was performed by 
one-step guanidinium isothiocyanate-acid phenol 
method (Chomezynski and Sacchi, 1987). The 
approach of reverse transcription (RT) followed 
by the polymerase chain reaction (PCR) (Rekik 
and Dea, 1994) was used to amplify the S1 region 
of the S genes, the entire HE genes, and entire 
ORF4 and ORF5 regions of the BCoV isolates 


studied. To access the error rate of the RT and 
Taq polymerase, clones from different RT-PCR 
events were sequenced. The primers were choosen 
according to the cDNA sequence of the Mebus 
strain of BCoV (GenBank accession numbers 
M31053, M30612 and M76372) (Abraham et al., 
1990a,b; Parker et al., 1990). The nt and deduced 
aa sequences were computer analysed with the aid 
of the GeneWorks 2,4 (IntelliGenetics Inc., 
Mountain View, Calif.) program. The nucleotide 
sequence accession numbers (EMBL/GenBank/ 
DDBJ) that have been assigned to the HE, Sl, 
ORF4 and ORF5 genomic regions of the EBCoV 
and RBCoV strains studied are: AF239307, and 
AF239316 for BCQ.1523; V06093 and AF239315 
for BCQ.571; AF239306, AHO010061 and 
AF239314 for BCQ.7373; L38962, AF239313 and 


AF239317 for BCQ.2590; AF239308 and 
AF239310 for BCO.43277; AF239309 and 
AF239311 for BCO.44175; AF339836 and 


AF239312 for BCQ.3994. To eliminate nt and aa 
changes that could be attributed to adaptation in 
HRT-G cells, comparison study was made with 
sequences of the HE and S genes of the avirulent 
Mebus strain initially propagated either in bovine 
fetal kidney cells or MDBK cells (GenBank acces- 
sion numbers M30612, D00662, S50936) (Abra- 
ham et al., 1990a,b; Kienzle et al., 1990). 
Comparison was also made with the sequences of 
the BCoV.Que strain, another attenuated NCD- 
BCoV strain isolated in 1972 from diarrheic calves 
in Quebec, Canada, which has been propagated 
for at least 41 passages in Vero, followed by 
10-15 passages in MDBK cells (Dea et al., 1980; 
Parker et al., 1990) prior sequencing studies 
(Rekik et al., 1994) (GenBank accession numbers 
D00662, M31053). 

The amplification of entire HE genes of RB- 
CoV strains OK-0514, BCO.43277, BCO.44175 
and BCQ.3994 were obtained using primers 
HE.Ad5 sens 5’ GGATCC GGATCC GCC GCC 
GCC ATG TTT TTG CTT CTT AGA TTT GTT 
C 3’ and HE.Ad5 antisens 5°GGATCC GGATCC 
TCA CTA AGC ATC ATG CAG CCT AGT 
ACC 3’, and comparison were made with the HE 
sequences of two Quebec WDBCoV strains 
(BCQ.2590 and BCQ.7373) and two Quebec 
NCDBCoV (BCQ.571 and BCQ.1523) previously 


Table 1 
The hemagglutinating (HA), esterase and RDE activities of selected Quebec field isolates of RBCoV and EBCoV, and their cross-reactivities in HAI tests to anti-HE 
MAbs directed against the WD-strain BCQ.2590* 


Strains Year isolation Disease HA titers? Esterase® activities (per 15 tl) RDE‘* HAI titers of anti-HE Mabs or rabbit antiserum® 
1D6-2 1D6-3.  9F2-IR —anti-Mebus 
BCoV.Meb 1972 NCD 2560 1.481 1 <20 <20 <20 1280 
BCQ.1523 1994 NCD 640 1.950 2 10 240 10240 1280 640 
BCQ.571 1989 NCD 640 1.782 1 <20 <20 <20 160 
BCQ.2590 1991 WD 2560 3.122 128 2560 2560 1280 160 
BCQ.7373 1992 WD 1280 2.181 256 10 240 5120 2560 640 
OK-0514-3 1996 Pn 2560 2.202 1 5120 5120 5120 320 
BCO.43277_ ~—- 1997 Pn 1280 2.178 1 5120 5120 40 320 
BCO.44175 1997 Pn 2560 2.170 1 5120 5120 80 320 
BCQ.3994 1998 Pn 5120 2.110 1 1280 1280 8 <20 320 


* BCoV, Bovine coronavirus; BCQ, Bovine coronavirus Quebec; BCO, Bovine coronavirus Ontario; Meb, reference avirulent Mebus strain; WD, Winter dysentery; 
NCD, Neonatal Calf Diarrhea; RBCoV, Respiratory bovine coronavirus; Pn, Pneumonitis. 

> Reciprocal of the highest dilution of infected cell culture fluids producing complete HA of rat erythrocytes after 1 h incubation at 4°C (averaged values for four 
tests after three to five passages on HRT-18 cells). 

© Absorbance at 405 nm after 5 min of reaction with 1 mM-p-nitrophenyl acetate: the extracellular virions of each BCoV isolates was purified from 100 ml of infected 
cell culture fluids following isopycnic ultracentrifugation on sucrose gradients. 

4 Reciprocal value of highest dilution producing complete disaggregation of BCoV-erythrocyte complexes after 1 h at 37°C. 

© Reciprocal values of highest MAb dilution inhibiting HA of rat erythrocytes by the various BCoV isolates (8 HA units) after 2 h at 4°C. 
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characterized in our laboratory (Dea et al., 1995; 
Kourtesis et al., 2001). Their nt and deduced aa 
sequences were then compared to those of the cell 
culture-attenuated reference © NCD-associated 
Mebus strain of BCoV which has been rese- 
quenced in our laboratory following passaged in 
HRT-G cells. For this reference strain, data were 
similar to those previously reported by other au- 
thors following propagation on the same cell line 
(Zhang et al., 1991b) or on FBK cells (Abraham 
et al., 1990a). Thus, HE gene sequences of the 
reference Mebus strain remained unchanged fol- 
lowing adaptation in HRT-G cells (Table 2). No 
nt substitution were also demonstrated when com- 
paring the HE genes of the attenuated Mebus and 
BCoV.Que strains. For all BCoV strains studied, 
the HE genes contained a large ORF of 1272 nts 
encoding a predicted protein of 424 amino acids 
(aa). Table 2 presents the aa changes observed for 
all strains analysed when compared to the refer- 
ence Mebus strain; the majority of the nucleotidic 
variations identified were silent mutations. The 
mutations of a leucine (L) by a proline (P) residue 
at aa position 5 and the substitution of a serine 
(S) by a proline (P) at aa position 367 determined 
for virulent BCoV isolates had been already re- 
ported (Zhang et al., 1991a; Dea et al., 1995). The 
mutation of an aspartic acid (D) by a glycine (G) 
at aa position 66 seems to be specific for all 
RBOCoV strains. 

For the S1 subunits of the S genes, a total of 
2304 nts were obtained for all BCoV strains stud- 
ied using primer pair S1.Ad5 sens 5‘'GGATCC 
GGATCC GCC GCC GCC ATG TTT TTG 
ATA CTT TTA ATT TCC 3’ and S1.Ad5 anti- 
sens AS 5'GGATCC GGATCC TCA TCT ACG 
ACT TCG TCT TTT TG 3’, corresponding to 
768 predicted aa residues. As previously reported 
by other investigators (Zhang et al., 1991a), se- 
quences of the Mebus and the BCoV.Que attenu- 
ated strains propagated either on FBK or MDBK 
cells were identical (data not shown). Most of the 
58 aa changes identified were strain-specific, but 
aa changes observed at position 11 (M changed 
for T) of the signal peptide and those observed at 
positions 40 (I for T), 169 (N for H), 173 (N for 
H), 248 (M for L), 253 (S for N), 465 (A for V), 
and 470 (D for H) could have resulted from the 


propagation on HRT-18 cells since they were also 
identified for the Mebus strain following propaga- 
tion on these cells (data not shown). The aa 
changes at positions 115 (D for K) and 173 (N for 
H) were observed for all virulent BCoV isolates. 
Eight aa_ substitutions, corresponding to aa 
residues 154, 182, 196, 215, 273, 320, 509 and 759, 
appeared to be specific to the EBCoV strains 
(Table 3). Only one of these variations, an as- 
paragine (N) residue replaced by a threonine (T) 
at aa position 509, was situated in the hypervari- 
able region of S1 (aa 452-593) (Rekik and Dea, 
1994). Another aa substitution at position 759, 
where an aspartic acid (D) was replaced by a 
valine (V) residue, appeared to be specific for 
virulent NCDBCoV strains. The inclusion of an 
asparagin (N) residue at aa position 320 specific 
to EBCoV strains did not result in the appearance 
of an additional putative glycosylation site, and 
the same situation was noted in the case of both 
WDBOoV strains for which a specific inclusion of 
an asparagin (N) residue at aa position 196 was 
revealed. 

Interestingly, eight aa substitutions in the S1 
subunit of the S gene also appeared to be more 
specific to RBCoV strains, corresponding to aa 
residues 118, 147, 256, 362, 510, 531, 578, and 691 
(Table 4). Three of these aa substitutions ap- 
peared to be situated in the hypervariable region 
of Sl. Three of these specific changes located at 
aa positions 510,531 and 578 had been already 
reported by other investigators (Chouljenko et al., 
1998) and were present in all RBCoV strains 
studied, except the Quebec RBCoV £ strain 
BCQ.3994 for which only the aa change at posi- 
tion 510 (S for T) was demonstrated. The aa 
substitution at positon 118 which was also previ- 
ously reported for the reference RBCoV strain 
OK-0514 was only demonstrated in the Ontario 
BCO.43277 strain. Three other aa changes iden- 
tified at aa positions 256 (M for L), 362 (S for T) 
and 691 (A for T) were observed only in the 
sequences of both Ontario RBCoV strains; the aa 
change at position 256 being also present in the 
Quebec RBCoV strain BCQ.3994. Chouljenko et 
al. (1998) had identified five other variations spe- 
cific to RBCoV strains situated at the N-terminus 
of the Sl portion of the S gene at aa positions 


Table 2 
Variations observed in the aa sequence of the HE glycoproteins of BCoV field strains compared to that of the reference Mebus strain* 


Amino acid Pathology and strains 


EBCoV strains RBCoV strains 
NCD WD 
Mebus Mebus BCoV.Que BCQ.571 BCQ.1523, BCQ.2590 BCQ.7373 OK-0514 BCO 43277 BCO 44175 BCQ 3994 
(HRT) (FBK) (MDBK) (HRT) (HRT) (HRT) (HRT) (ART) (HRT) (HRT) (HRT) 
5 L x sf P P P Pp Pp p Pp i) 
49 N * * 7 T = T T T T T 
66 D bs * Es * * * G G G G 
245 F * * * * * L * * * * 
282 D * * * * * G * * * * 
344 D * * * * A % * * * * 
367 S * * P P P P P P P P 
392 L 7 * I I * I I I I I 


@ An asterisk (*) indicates a consensus with the reference Mebus strain. 
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Table 3 
Variations observed in the aa sequence of the S1 portion of the S glycoproteins of BCoV field strains associated to WD, NCD, or pneumonia, specific to virulent 
EBCoV strains when compared to the reference Mebus strain* 


Amino acid Pathology and strains 


EBCoV strains RBOCoV strains 
NCD WD 
Mebus Mebus BCoV.Que BCQ.571 BCQ.1523, = BCQ.7373 ~=BCQ.2590 OK-0514 BCO 43277 BCO 44175 BCQ 3994 
(HRT) (FBK) (MDBK) (ART) (ART) (HRT) (ART) (HRT) (ART) (HRT) (HRT) 
182 E * * K K K K * * * * 
215 E * . K K kK K * * * * 
273 A * P P P P a * * 
320 D * * N N N N - * * * 


* An asterisk (*) indicates a consensus with the reference Mebus strain. 
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Table 4 
Variations observed in the aa sequence of the S1 portion of the S glycoproteins of BCoV field strains associated to WD, NCD, or pneumonia, specific to virulent 
RBCoV strains when compared to the reference Mebus strain* 


Amino acid Pathology and strains 


EBCoV strains RBCoV strains 
NCD WD 
Mebus Mebus BCoV.Que BCQ.571 BCQ.1523, = BCQ.2590 BCQ.7373 OK-0514 BCO 43277 BCO 44175 BCQ 3994 
(HRT) (FBK) (MDBK) (HRT) (HRT) (HRT) (HRT) (HRT) (HRT) (HRT) (HRT) 
118 M * * ba * K * K K M M 
510 S is * ba ig = = T T T T 
531 N = * ia G * = G G G ig 


* An asterisk (*) indicates a consensus with the reference Mebus strain. 


LS-E€F (1007) 9L Yyotvasay sn.it4 / ‘JD Ja spUaD "W-V 


Is 


By? A.-M. Gélinas et al. / Virus Research 76 (2001) 43-57 


455) (05/9 6) 901011 


Lom <— 1014 bp 


—_—— — oe oe 
os. + 638 bp 


Fig. 1. RT-PCR amplification of the genomic region between 
the spike (S) and matrix (M) protein genes of BCoV isolates 
associated to different clinical syndromes as compared to the 
corresponding region of the porcine HEV and human HCoV- 
OC43 viruses. As expected, the use of primer pairs 
BCoV.5327S and BCoV.ORF4AS yielded DNA amplicons of 
1014 bp for BCoV isolates associated to either enteric or 
respiratory diseases, whereas smaller DNA amplicons of 638 
bp were obtained for both HEV-67N and HCoV-OC43 
viruses. L, 1 kb DNA ladder; 1, Mebus strain (NCD); 2, BCQ. 
2590 (WD); 3, BCQ.7373 (WD); 4, BCQ.571 (NCD); 5, 
BCQ.1523 (NCD); 6, HEV-67N; 7, HCoV-OC43; 8, OK-0514 
(RBCoV); 9, BCO.43277 (RBCoV); 10, BCO.44175 (RBCoV); 
11, BCQ.3994 (RBCoV). 


11,115, 118,173 and 179. In the present study, 
these subtitutions were also present in the Sl 
sequence of EBCoV strains. 

The genomic region between the S and M 
genes, corresponding to the ORF4 and ORF5 
genes, was amplified by RT-PCR from NCD- and 
WD.-associated BCoV, as well as from the four 
RBCoV strains studied, using primer pair 
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BCoV.5327 sens 5'GGATCC GGATCC ATG 
TGG TGG TTG TTG TGA TGA 3’ and 
BCoV.ORF4 antisens 5°'GATTC GAATTC AAC 
GTGC ATC CAC ATC AAG AAC 3’. For all 
BCoV isolates studied, DNA amplicons of 1014 
base pairs were obtained, whereas smaller DNA 
amplicons of 638 base pairs were obtained for 
both HCoV-OC43 and HEV-67N strains (Fig. 1). 
For both types of EBCoV strains, the ORF4 
region apparently encodes for two small non- 
structural proteins of 4.9 and 4.8 kDa and com- 
prises a total of 302 nts. The alignments of the aa 
sequences of the ORF4 region encoding the 4.9 
kDa protein of all the BCoV strains studied are 
presented in Fig. 2. These sequences were com- 
pared to that of the reference Mebus strain and to 
those of HEV-NT9, HEV 67N and HCoV-OC43 
strains (GenBank accession numbers X89861, 
X89863 and M99576). As previously reported by 
other authors (Mounir and Talbot, 1993; Vieler et 
al., 1996), the HEV and HCoV-OC43 strains 
showed a major deletion in the genomic region 
corresponding to the ORF4, which was not the 
case for EBCoV and BCoV strains. However, a 
specific non-sense codon was identified at nt posi- 
tion 88 of the ORF4 of the four RBCoV isolates 
studied resulting in a premature termination of 
the protein which comprises 29 rather than 43 aa 
residues. In contrast to data obtained by previous 
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Fig. 2. Comparison of aa sequences of the ORF4 region encoding the 4.9 kDa protein. In comparison to BCoV enteropathogenic 
strains, the human coronavirus OC43 and the porcine HEV ORF 4 lack two open reading frames, whereas all RBCoV field isolates 
tested possess these two open reading frames but the 4.9 ns protein is apparently not expressed because of the introduction of an 
antisense codon at position 88. Sequences of strains HEV-NT9, HEV-67N and HCoV-OC43 correspond to GenBank accession 


numbers X89861, X89863, and M99576, respectively. 
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investigators (Chouljenko et al., 1998), all virulent 
field strains, and not only RBCoV strains, con- 
tained a deletion (TAAT) located 31-34 bp up- 
stream from the AUG of the 12.7 kDa ORF in 
the genomic region between the 4.8 and 12.7 kDA 
ORFs (data not shown). Only a four strain-spe- 
cific aa changes were observed in the ORFS re- 
gion encoding the 12.7 and 9.5 kDa (E) proteins. 
The only substitution determined within the E 
gene was located at aa position 53 where a valine 
was replaced by a glycine for all virulent strains, 
as previously described (Mounir and Talbot, 
1993; Chouljenko et al., 1998). 

The data obtained in the present study indicate 
that genomic, biologic and antigenic variations, as 
detected by MAbs directed to epitopes of two 
different antigenic domains of the HE glyco- 
protein of a reference WDBCoV strain (Milane et 
al., 1997) exist among EBCoV and RBCoV 
strains. However, the majority of aa changes de- 
tected seems to be scattered over different struc- 
tural genes and appeared to be strain-specific, in 
agreement with previous findings on NCD- and 
WD.-associated strains (Zhang et al., 1991a,b; 
Rekik and Dea, 1994; Chouljenko et al., 1998; 
Kourtesis et al., 2001). Sequencing analyses were 
restricted to genes encoding the envelope-associ- 
ated HE, S (S1 portion) and E proteins of the 
BCoV virion, as well as putative ns proteins 
(4.9, 4.8 and 12.7 kDa) encoded between the spike 
and membrane protein genes, which have been 
previously incriminated in the tissue tropism 
(Mounir and Talbot, 1993; Vieler et al., 1996), 
pathogenesis, humoral and cellular immune re- 
sponses (Spaan et al., 1988 (a review)), as well as 
to genomic and antigenic variabilities amongst 
BCoV field isolates (El Ghorr et al., 1989; 
Michaud and Dea, 1993; Rekik and Dea, 1994; 
Tsunemitsu and Saif, 1995; Milane et al., 1997; 
Chouljenko et al., 1998; Hasoksuz et al., 1999). 
The N and M protein genes were not analysed 
since previously serological and genomic compari- 
son studies revealed that the latter two structural 
proteins are well conserved among field strains of 
BCoV associated to various clinical conditions, 
and amongst members of the hemagglutinating 
coronaviruses subgroup (Lapps et al., 1987; 
Schockley et al., 1987; Verbeek and Tijssen, 1991; 


Mounir and Talbot, 1992; Chouljenko et al., 
1998). A recent study demonstrated that se- 
quences of the N were more conserved among 
RBCoV strains than to EBCoV strains compared, 
although the few aa changes identified were 
strain-specific (Chouljenko et al., 1998). 

Data further demonstrated that both envelope- 
associated glycoproteins, particularly the S1 sub- 
unit of the S peplomeric glycoprotein, appeared to 
contain aa changes specific to either RBCoV or 
EBCoV strains, in agreement with recent findings 
on RBCoV strains characterized in Louisiana 
(Chouljenko et al., 1998). Because the S protein is 
more efficient that the HE protein in hemaggluti- 
nation (Vlasak et al., 1988; Schultze et al., 1991b), 
it is expected that antigenic variability would have 
been more evident using anti-S MAbs, rather than 
anti-HE MAbs. Interestingly, none of the anti-HE 
MAbs showed reactivity by HAI and indirect 
immunofluorescence (data not shown) to the 
reference Mebus strain, as well as two other 
NCDBCoV isolates, as previously reported (Mi- 
lane et al., 1997). However, these anti-HE MAbs 
could not be used for specific identification of 
WDBCoV strains since they all reacted to the 
RBCoV field isolates. On the other hand, all 
BCoV strains tested reacted positively by HAI 
against both rabbit hyperimmune sera (anti- 
Mebus and anti-WD-2590). This may be ex- 
plained by the fact that hemagglutination through 
S may still occur when the HE hemagglutinin is 
inhibited. However, epitopes of the S associated 
to virus neutralization are not necessarily associ- 
ated to those associated to HAI, since both rabbit 
hyperimmune antisera were devoid of neutraliza- 
tion activities, in agreement with previous findings 
(Vautherot et al., 1992; Michaud and Dea, 1993). 

One important finding which has been also 
observed previously for the reference WDBCoV 
strain BCQ.2590 (Milane et al., 1997) is the 
highest enzymatic activities (RDE and esterase 
activities) of WDBCoV field isolates in compari- 
son to that of NCDBCoV and RBCoV isolates. 
Although earlier studies have demonstrated that 
BCoV esterase activity is probably not directly 
implicated in the processus of infection, particu- 
larly virus uptake (Vlasak et al., 1988; Schultze et 
al., 1991a), it could be related to the virulence by 
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facilating elution of the mature virions from the 
surface of infected cells, thus a possible higher 
potential for rapid dissemination (higher conta- 
giosity) in the herds. Results from comparative 
studies of the HA properties of NCD- and WD- 
associated BCoV isolates also suggest that the 
affinity for N-acetyl-9-O neuraminic acid-contain- 
ing receptors is probably dependent on the RDE 
activity of the HE protein (Dea et al., 1995). 
Accordingly, outbreaks of WD are usually of 
short duration and affected rapidly most of the 
animals within the herd, in contrast to NCD 
outbreaks that rather tend to evolve to a chronic 
status since newborn calves can shed the virus in 
their feces for long periods after recovery (Crouch 
et al., 1985; Dea et al., 1995; Fukutomi et al., 
1999). It has been previously suggested that en- 
hancement of the RDE activity of WDBCoV 
strains could be in part attributed to conforma- 
tional changes of the HE protein due to inclusion 
of proline residues in the signal peptide and in the 
vicinity of the sequences of the putative AE do- 
main (FGDS) of the HE protein. Data of the 
present study do not support this hypothesis since 
both proline substitutions were identified for all 
virulent BCoV strains studied, including those 
characterized by other investigators (Zhang et al., 
1991b). The substitution of an aspartic acid (D) 
by a glycine (G) at aa position 66 of the HE 
glycoprotein was also previously identified as a 
specific marker for RBCoV strains (Chouljenko et 
al., 1998). 

As also previously reported (Chouljenko et al., 
1998), the region located between the S and M 
protein genes appeared of a particular interest 
when comparing EBCoV and RBCovV isolates. 
Indeed, this region contained four potential ORFs 
for all RBCoV isolates studied; thus BCoV iso- 
lates associated to different clinical symptoms did 
not have the 302 nt deletion previously identified 
for the human coronavirus HCoV-OC43 strain 
and HEV strains of swine that lack the two ORFs 
which potentially encode two non-structural (ns) 
proteins of 4.9 and 4.8 kDa (Mounir and Talbot, 
1993; Vieler et al., 1996). Their absence in cases of 
both viruses indicates that they are most probably 
not essential for virus replication, at least in 
HRT-18 cells. Interestingly, an antisense mutation 


resulting in the truncation of the ns 4.9 protein 
was identified only for RBCoV field isolates, 
which may indicate that this protein could be 
involved, as for HEV and HCoV-OC43, in the 
preferential respiratory tropism. Further investi- 
gations are required to verify this hypothesis. This 
situation contrasts considerably with that of the 
porcine transmissible gastroenteritis virus (TGEV) 
for which deletions within the first 250 aa of S1 
portion of S are associated with respiratory 
tropism and reduced enteropathogenicity (Wesley 
et al., 1991). On the other hand, the genomic 
regions encoding for the putative 12.7 kDa ns 
protein and the 9.5 kDa membranous protein E 
are probably essential for virus replication and 
infection process since both proteins appeared 
highly conserved among BCoV strains associated 
to various clinical symptoms. 

Although differences were found by sequencing 
analysis and patterns could be correlated some- 
what with pathogenic properties of the BCoV 
strains characterized, it remains that the region of 
the genome studied represents only 25% of the 
entire BCoV genome. Therefore, data obtained do 
not eliminate the fact that differences among 
strains that might give rise to differences in patho- 
genic properties could lie anywhere within the 
genome, including non-coding regions. For exam- 
ple, previous investigators have demonstrated the 
existence of nt variations in the genomic region 
located in the N-terminus of the ORF2 encoding 
a 32 kDa non-structural protein, highly expressed 
in BCoV-infected cells, or in cells infected by 
other members of the hemagglutinating coro- 
naviruses (Cox et al., 1989, 1991; Labonté et al., 
1995). In comparison to the EBCoV attenuated 
BCoV.Que strain, two virulent RBCoV strains 
isolated in Louisiana, as well as the virulent LY- 
138 EBCoV strain, have been found to contain 
two frame-shift mutations which resuted in a nine 
aa segment near the carboxy-terminus which was 
different from the corresponding aa sequence spe- 
cified by the avirulent BCoV.Que - strain 
(Chouljenko et al., 1998). Sequencing of the 32 
kDa non-structural protein of the attenuated and 
virulent Mebus stains, as well as the same region 
of other virulent strains should substantiate these 
differences between virulent and avirulent strains. 
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Other regions may also be involved, namely the 
large ORF1 region encoding the RNA-poly- 
merase, replicase, helicase and other proteases 
involved in virus replication and maturation 
(Spaan et al., 1988). Such findings have been 
recently demonstrated in the case of the porcine 
reproductive and respiratory syndrome virus, a 
member of the Arteriviridae family recently 
grouped in the order Nidovirales, together with 
members of the Coronaviridae family (De Vries et 
al., 1997). The genome length of Arterivirus is 15 
kb rather than 28 kb, which makes genomic com- 
parison studies more easier than in case of coro- 
naviruses, particularly for analysis of the ORF1 
gene representing near 66% of the entire gemone 
for both types of viruses. The ORF1 is separated 
into ORFla and ORF1b regions by a ribosomal 
frameshift; six putative Nsps (non-structural 
proteins) are predicted to result from ORFla 
polyprotein processing, whereas ORF 1b is cleaved 
by ORFla proteases into four proteins: RdPp 
(RNA dependent RNA polymerase) and putative 
cleavage products CP2, CP3 and CP4 (De Vries et 
al., 1997; Allende et al., 1999a). The Nsp2 protein 
has been shown highly variable among. ar- 
teriviruses, e.g. PRRSV (740 and 861 residues 
long for European and North American strains, 
respectively), LDV (733 residues long) and EAV 
(401 residues long), with similarity only in the 
amino- and carboxy-terminal domains, and _ is 
thought to be involved in species-specific func- 
tions. The low degree of identity and the size 
differences between Nsp2 of European and North 
American isolates of PRRSV suggests that vari- 
ability within this protein may have other than 
species-specific functions. Whether this involves 
biological differences observed between virus iso- 
lates, such as the degree of pathogenicity in the 
swine host, remains to be determined. Further- 
more, aa changes distributed in the different re- 
gions of the ORF1 gene have been identified 
which may provide the molecular bases for loss of 
virulence of the live-attenuated PRRSV vaccines 
(Allende et al., 1999b). To our knowledge, such 
variability in the various regions of the ORF1 of 
BCoVs associated with different clinical symp- 
toms has not been yet demonstrated, but could 
give rise to differences in pathogenic properties, 
tropism and antigenic variabilities of field isolates. 


Acknowledgements 


This research was supported by an operating 
grant received from the National Sciences and 
Engineering Research Council of Canada (Grant 
# OGPOO46412). The authors thank Nicole 
Sawyer and Louise Wilson for technical assis- 
tance. The authors greatly acknowledge J. Storz 
(Veterinary School of Medecine, Louisiana’s State 
University) for providing HRT-18G cells and 
OK-0514 RBCoV reference strain, and S. Carman 
(Ontario Ministry of Agriculture and Fisheries, 
Ontario Veterinary College, Guelph, Ontario) for 
providing both ontario RBCoV strains. We are 
also grateful to Dr C. Montpetit, Direction de la 
Santé Animale, Ministére de lAgriculture, des 
Pécheries et de l’Alimentation du Québec, for 
providing field strains of bovine coronavirus. 


References 


Abraham, S., Kienzle, T.E., Lapps, W., Brian, D.A., 1990a. 
Deduced sequence of the bovine coronavirus spike protein 
and identification of the internal proteolytic cleavage site. 
Virology 176, 296-301. 

Abraham, S., Kienzle, T.E., Lapps, W., Brian, D.A., 1990b. 
Sequence and expression analysis of potential nonstruc- 
tural proteins of 4.9, 4.8, 12.7 and 9.5 kDa encoded 
between the spike and membrane protien genes of the 
bovine coronavirus. Virology 177, 488—495. 

Allende, R., Lewis, T.L., Lu, Z., Rock, D.L., Kutish, G.F., 
Ali, A., Doster, A.R., Osario, F.A., 1999a. North Ameri- 
can and European porcine reproductive and respiratory 
syndrome viruses differ in non-structural protein coding 
regions. J. Gen. Virol. 80, 307-315. 

Allende, R., Kutish, G.F., Laegreid, W., Lu, Z., Lewis, T.L., 
Rock, D.L., Friesen, G.F., Galeota, J.A., Dosyer, A.R., 
Osario, F.A., 1999b. Mutations in the genome of porcine 
reproductive and respiratory syndrome virus responsible 
for the attenuation phenotype. Arch. Virol. 145, 1149— 
1161. 

Benfield, D.A., Saif, L.J., 1990. Cell culture propagation of a 
coronavirus isolated from cows with winter dysentery. J. 
Clin. Microbiol. 28, 1454-1457. 

Chomezynski, P., Sacchi, N., 1987. Single-step method of 
RNA isolation by acid guanidium thiocyanate—phenol— 
chloroform extraction. Ann. Biochem. 162, 156-159. 

Chouljenko, V.N., Kousoulas, K.G., Lin, X., Storz, J., 1998. 
Nucleotide and predicted amino acid sequences of all genes 
encoded by the 3’ genomic portion (9.5 kb) of respiratory 
bovine coronaviruses and comparisons among respiratory 
and enteric coronaviruses. Virus Genes 17, 33-42. 


56 A.-M. Gélinas et al. / Virus Research 76 (2001) 43-57 


Cox, G.L., Parker, M.D., Babiuk, L.A., 1989. The sequence of 
cDNA of bovine coronavirus 32K nonstructural gene. 
Nucleic Acids Res. 17, 5847. 

Cox, G.L., Parker, M.D., Babiuk, L.A., 1991. Bovine coro- 
navirus nonstructural protein ns2 is a phosphoprotein. 
Virology 196, 509-512. 

Crouch, C.F., Bielefeldt, O.H., Watts, T.C., Babiuk, L.A., 
1985. Chronic shedding of bovine enteric coronavirus anti- 
gen—antibody complexes by clinically normal cows. J. Gen. 
Virol. 66, 1489-1500. 

Dea, S., Roy, R.S., Begin, M.E., 1980. Bovine coronavirus 
isolation in continuous cell lines. Am. J. Vet. Res. 41, 
30-38. 

Dea, S., Garzon, S., Strykowski, H., Tijssen, P., 1989. Ultra- 
structure and protein A-gold immunolabelling of HRT-18 
cells infected with turkey enteric coronavirus. Vet. Micro- 
biol. 20, 21-33. 

Dea, S., Verbeek, A.J., Tijssen, P., 1990. Antigenic and ge- 
nomic relationships among turkey and bovine enteric coro- 
naviruses. J. Virol. 64, 3112-3118. 

Dea, S., Michaud, L., Milane, G., 1995. Comparison of bovine 
coronavirus isolates associated with neonatal calf di- 
arrhoea and winter dysentery in adult dairy cattle in Qué- 
bec. J. Gen. Virol. 76, 1263-1270. 

De Vries, A.A.F., Horzinek, M.C., Rottier, P.J.M., De Groot, 
R.J., 1997. The genome organization of the Nidovirales: 
similarities and differences between Arteri-, Toro-, and 
Coronaviruses. Sem. Virol. 8, 33—47. 

El Ghorr, A.A., Snodgrass, D.R., Scott, F.M.M., Campbell, 
L, 1989. A serological comparison of bovine coronavirus 
strains. Arch. Virol. 104, 241-248. 

Fukutomi, T., Tsunemitsu, H., Akashi, H., 1999. Detection of 
bovine coronaviruses from adult cows with epizootic di- 
arrhea and their antigenic and biological diversities. Arch. 
Virol. 144, 997-1006. 

Hasoksuz, M., Lathrop, S., Al-dubaib, M.A., Lewis, P., Saif, 
L.J., 1999. Antigenic variation among bovine enteric coro- 
naviruses (BECV) and bovine respiratory coronaviruses 
(BRCV) detected using monoclonal antibodies. Arch. Vi- 
rol. 144, 2441-2447. 

Hogue, B.G., Brian, D.A., 1986. Structural proteins of human 
respiratory coronavirus OC43. Virus Res. 5, 131-144. 
Kienzle, T.E., Abraham, S., Hogue, B.G., Brian, D.A., 1990. 
Structure and expression of the bovine coronavirus hemag- 

glutinin protein. Adv. Exp. Mol. Biol. 276, 95-102. 

King, B., Poots, B.J., Brian, D.A., 1985. Bovine coronavirus 
hemagglutinin protein. Virus Res. 2, 53-59. 

Kourtesis, A.B., Gelinas, A.M., Dea, S., 2001. Genomic and 
antigenic variations of the HE glycoprotein of bovine 
coronaviruses associated with neonatal calf diarrhea and 
winter dysentery. Arch. Virol. (in press). 

Labonté, P., Mounir, S., Talbot, P., 1995. Sequence and 
expression of the ns2 protein gene of human coronavirus 
0C43. J. Gen. Virol. 76, 431-435. 

Lai, M.M.C., 1986. Coronavirus leader-RNA-primed tran- 
scription: an alternative mechanism to RNA splicing. 
Bioassay 5, 257-260. 


Lapps, W., Hogue, B.G., Brian, D.A., 1987. Sequence analysis 
of the bovine coronavirus nucleo-capsid aand matrix 
protein. Virology 157, 47—57. 

Mebus, C.A., Stair, E.L., Rhodes, M.B., Twiehaus, M.J., 
1973. Neonatal calf diarrhea: propagation, attenuation and 
characteristics of a coronavirus-like agent. Am. J. Vet. Res. 
34, 145-150. 

Michaud, L., Dea, S., 1993. Characterization of monoclonal 
antibodies to bovine enteric coronavirus and antigenic 
variability among Quebec isolates. Arch. Virol. 131, 455— 
465. 

Milane, G., Kourtesis, A.B., Dea, S., 1997. Characterization 
of monoclonal antibodies to the hemagglutinin—esterase 
glycoprotein of a bovine coronavirus associated with win- 
ter dysentery and cross-reactivity to field isolates. J. Clin. 
Microbiol. 35, 33—40. 

Mounir, S., Talbot, P.J., 1993. Human coronavirus OC43 
RNA 4 lacks two open reading frames located downstream 
of the S gene of bovine coronavirus. Virology 192, 355— 
360. 

Mounir, S., Talbot, P.J., 1992. Sequence analysis of the mem- 
brane protein gene of human coronavirus OC43 and evi- 
dence for O-glycosylation. J. Gen. Virol. 73, 2731-2736. 

Parker, M.D., Yoo, D., Babiuk, L.A., 1990. Expression and 
secretion of the bovine coronavirus hemagglutinin esterase 
glycoprotein by insect cells infected with recombinant bac- 
uloviruses. J. Virol. 64, 1625-1629. 

Rekik, M.R., Dea, S., 1994. Comparative sequence analysis of 
a polymorphic region of the spike glycoprotein S1 subunit 
of enteric bovine coronavirus isolates. Arch. Virol. 135, 
319-331. 

Reynolds, D.J., Debney, T.J., Hall, G.A., Thomas, L.H., 
Parsons, K.R., 1985. Studies on the relationship between 
coronaviruses from the intestinal and respiratory tracts of 
calves. Arch. Virol. 85, 71-83. 

Schockley, L.J., Kapke, P.A., Lapps, W., Brian, D.A., Potgi- 
eter, L.N.D., Woods, R., 1987. Diagnosis of porcine and 
bovine enteric coronavirus infections using cloned cDNA 
probes. J. Clin. Microbiol. 25, 1591-1596. 

Schultze, B., Wahn, K., Klenk, H.D., Herrler, G., 1991a. 
Isolated HE-protein from hemagglutinating — en- 
cephalomyelitis virus and bovine coronavirus has receptor- 
binding activity. Virology 180, 221-228. 

Schultze, B., Gross, H.J., Brossmer, R., Herrler, G., 1991b. 
The S protein of bovine coronavirus is a hemagglutinin 
recognizing 9-O-acetylated sialic acid as a receptor deter- 
minant. J. Virol. 65, 6232—6237. 

Spaan, W.J.M., Rottier, P.J.M., Horzinek, M.C., Van Der 
Zeijst, B.A.M., 1981. Isolation and identification of virus- 
specific mRNAs in cells infected with mouse hepatitis virus 
(MHV-AS59). Virology 108, 581-592. 

Spaan, W., Cavanagh, D., Horzinek, M.C., 1988. Coro- 
naviruses: structure and genome expression. J. Gen. Virol. 
69, 2939-2952. 

Storz, J., Zhang, X.M., Rott, R., 1992. Comparison of hemag- 
glutinating, receptor destroying, and acetylesterase activi- 
ties of avirulent and virulent bovine coronavirus strains. 
Arch. Virol. 125, 193-204. 


A.-M. Gélinas et al. / Virus Research 76 (2001) 43-57 57 


Storz, J., Stine, L., Liem, A., Anderson, G.A., 1996. Coro- 
navirus isolation from nasal swabs samples in cattle with 
signs of respiratory tract disease after shipping. J. Am. Vet. 
Med. Assoc. 208, 1452-1454. 

Storz, J., Lin, X., Purdy, C.W., Chouljenko, V.N., Kousoulas, 
K.G., Enright, F.M., Gilmore, W.C., Briggs, R.E., Loan, 
R.W., 2000. Coronavirus and Pasteurella infections in 
bovine shipping fever pneumonia and Evans’ criteria for 
causation. J. Clin. Microbiol. 38, 3291-3298. 

Tsunemitsu, H., Saif, L.J., 1995. Antigenic and biological 
comparisons of bovine coronaviruses derived from neona- 
tal calf diarrhea and winter dysentery of adult cattle. Arch. 
Virol. 140, 1303-1311. 

Vautherot, J.F., Madelaine, M.F., Boireau, P., Laporte, J., 
1992. Bovine coronavirus peplomer glycoproteins: detailed 
antigenic analysis of S1, S2 and HE. J. Gen. Virol. 73, 
1725-1737. 

Verbeek, A., Tijssen, P., 1991. Sequence analysis of the turkey 
enteric coronavirus nucleocapsid and membrane protein 
genes: a close genomic relationship with bovine coro- 
navirus. J. Gen. Virol. 72, 1659-1666. 


Vieler, E., Schlapp, T., Herbst, W., 1996. The region between 
the M and S genes of porcine hemagglutinating en- 
cephalomyelitis virus is highly similar to human coro- 
navirus OC43, J. Gen. Virol. 77, 1443-1447. 

Vlasak, R., Luytjes, W., Leider, J., Spaan, W., Palese, P., 
1988. The E3 protein of bovine coronavirus is a receptor- 
destroying enzyme with acetylesterase activity. J. Virol. 62, 
4686-4690. 

Wesley, R.D., Woods, R.D., Cheung, A.K., 1991. Genetic 
analysis of porcine respiratory coronavirus, an attenuated 
variant of transmissible gastroenteritis virus. J. Virol. 65, 
3369-3373. 

Zhang, X., Kousoulas, K.G., Storz, J., 1991a. Comparisson of 
the nucleotide and deduced amino acid sequences of the S 
genes specified by virulent strains of bovine coronaviruses. 
Virology 183, 397-404. 

Zhang, X., Kousoulas, K.G., Storz, J., 1991b. The hemagglu- 
tinin/esterase glycoprotein of bovine coronaviruses: se- 
quence and functional comparisons between virulent and 
avirulent strains. Virology 185, 847-852. 


